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ABSTRACT 


1 

This  report  deals  with  the  use  of  high  frequency  asymptotics 
to  evaluate  and  possibly  facilitate,  the  convergence  of  the 
transition  matrix  formulation  of  scattering.  For  illustrative 
purposes,  elastic  and  acoustic  scattering  from  spherical  and 
spheroidal  cavities  is  considered  in  some  detail.  References 
are  presented  for  other  types  of  inclusions.  Finally,  a  general 
method  is  suggested  for  integrating  high  frequency  asymptotic 
solutions  into  the  T-matrix  formulation  itself  with  the  purpose 
of  improving  convergence.  The  method  is  implemented  for  the 
case  of  a  sphericaL,  and  a  prolate  spheroidal,  cavity. 
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I.  INTRODUCTION 

The  transition  matrix  (T-matrix)  formulation  of  scattering  is  the 

most  recent  addition  to  the  well-known  classical  methods:  separation  of 

1  2 

variables,  the  integral  equation  method,  and  the  variational  method.  ' 

This  method  has  been  applied  to  acoustic,  elastic,  and  electromagnetic 
media  and  to  a  variety  of  complicated  geometries. 

All  of  the  above  methods,  including  the  T-matrix  method,  have  the 
drawback  of  poor  convergence  when  the  ratio  D/X  >>  1,  where  D  is  the 
characteristic  dimension  of  the  scatterer  and  X  is  the  characteristic 
wavelength  of  the  scattered  field.  Typically,  the  number  of  terms 
required  for  convergence  is  of  order  2D/X.  However,  "coordinate  stretching" 
can  be  used  to  reduce  the  extent  of  the  calculations  implied  by  this 

7 

criterion.  While  the  TM  method  appears  to  be  a  most  efficient  formulation, 
and  consequently,  the  most  widely  studied  over  these  last  few  years,  it 
too  entails  numerical  difficulties  at  high  frequencies. 

For  problems  in  steady  state  elastic  and  acoustic  scattering,  the 

character  and  functional  dependence  of  the  scattered  field,  as  well  as 

its  amplitude,  vary  with  frequency.  In  the  low  frequency  or  Rayleigh 

region,  defined  by  D/X  <  1,  the  scattered  field  is  completely  defined  in 

terms  of  only  the  gross  characteristics  of  the  scatterer,  viz.  it’s  overall 

8 

dimensions,  relative  density  and  compressibility.  Computationally,  this 
region  generally  involves  little  difficulty  using  any  of  the  traditional 
analytical  techniques.  The  low  frequency  regime  is  followed  by  the 
resonance,  or  standing  wave,  region  whereby  the  scatterer  dimensions  are 
comparable  to  the  characteristic  wavelength (s) ,  i.e.,  D/X  =  0(1-10).  It 
is  at  the  upper  end  of  this  regime  where  numerical  formulations,  either 
modal  or  integral,  begin  to  exhibit  some  computational  difficulty;  lack 
of  convergence,  round  off  errors,  the  existence  of  forbidden  frequencies, 
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As  frequency  increases  further,  or  correspondingly  as  wavelength 
decreases,  the  scatterer  begins  to  measure  many  wavelengths.  Here  the 
scattered  field  is  more  conveniently  described  in  terms  of  travelling 
waves  which  are  generated  at  the  boundary  of,  and  within,  the  scatterer. 
These  travelling  waves,  because  of  diffraction  and  dissipation  are 
attenuated  as  they  circumnavigate  the  scatterer.  For  present  purposes 
this  is  defined  as  the  creeping  wave  region.*  Finally,  the  limiting  high 
frequency  asynqptotic  (Kirchhoff  or  Physical  Optics)  reqion  is  defined  as 
that  region  where  ray  theory  applies,  i.e.,  where  the  scatterer  can  be 
considered  to  be  locally  reacting  and  where  each  scattering  element  acts 
as  if  it  were  in  an  infinite  baffle. 

It  is  in  these  latter  two  frequency  regimes  where  the  computational 
difficulties  encountered  by  numerical  formulations,  including  the  T  matrix, 
become  severe.  The  purpose  of  this  study  is  to  analyze  the  potential  role 
that  high  frequency  asymptotic  scattering  formulations  may  play  in  serving 
(a)  as  a  convergence  check  for  T  matrix  solutions,  and  (B)  as  a  means  for 
improving  the  convergence  of  T  matrix  formulations.  For  discussion 
purposes  only  elastic  and  acoustic  media,  and  spherical  or  spheroidal 
scatterers,  are  considered.  Also,  the  high  frequency  formulations  are 
limited  to  the  Kirchhoff  and  creeping  wave  type. 

In  Section  II  the  problem  of  elastic  scattering  from  spherical 
inclusions  is  reviewed,  together  with  the  corresponding  T-matrix  and 
wave  harmonic  formulations.  The  corresponding  high  frequency  asymptotic 
solutions  are  also  discussed  in  Section  II.  In  Section  III  a  general 
formalism  for  using  these  high  frequency  solutions  to  complement  the 
T-matrix  formulation  is  presented  together  with  illustrative  examples. 
Specifically,  the  examples  are  scalar  scattering  from  a  spherical,  and  a 
prolate  spheroidal  cavity. 


*  It  is  noted  that  this  definition  of  the  creeping  wave  region  is  more 
restrictive  than  need  be  since,  in  theory,  creeping  wave  solutions  can 
be  as  general  as  is  the  corresponding  wave  harmonic  series  from  which 
they  are  derived. 
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II.  ELASTIC  SCATTERING  PROM  SPHERICAL  INCLUSIONS 
A.  Governing  Equation*  of  Motion 

The  harmonic  response  of  a  three  dimensional  linearly  elastic 
solid  can  be  fully  defined  in  terms  of  a  scalar  and  a  vector  potential 
which  satisfy  a  scalar  and  vector  wave  equation  respectively.  In 
spherical  coordinates  (Pig.  1),  the  vector  wave  equation  is  itself 
separable  and  one  is  left  with  the  following  governing  equations  of 
action 


2  2 

c^V  4  u>  ♦  “  0 

2  2  2 

c^v  +  w  t|i  -  o 


(la) 

(lb) 


2  2  2 
c#V  x  +  u  x  ■  o 


(lc) 


from  which  the  associated  displacement  field  may  be  calculated.  Using 
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spherical  coordinates 


u  ■  ii  +  a  ■■■  -  rV2v 


3r 


(2a) 


3r 


_  1 

U6  “  r  30  rsin0 


1  3(r»)  1  3  (rx) 

30  r  303r 


(2b) 


1  ii  .  I  3  (rid)  +  a  1 

rsinO  30  r  36  rsinO 


303r 


(2c) 


In  the  above  equations  u>  is  the  circular  frequency  and  c^  and  cg  represent 
the  dilatational  and  shear  wave  speeds,  respectively.  The  corresponding 

9 

stress  field  can  also  be  explicitly  determined.  For  an  axisynsietric, 
i.e.,  3  independent,  problem  such  as  is  considered  in  this  report,  ■>  0. 

In  the  case  of  a  cavity  of  radius  a  centered  at  the  origin  the  boundary 
conditions  to  be  satisfied  are 

a  (a, 0,6)  -  a  .(a, 0,6)  -  a  .(a, 0,8)  -  0  (3a) 

rr  ro  rs 

For  a  rigid  and  imovable  inclusion,  also  of  radius  a  and  centered  at  the 
origin,  the  boundary  conditions  are 
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ur<a,0,8)  -  u0(a,8,g>  -  0  (3b) 

A  dilatational  source  may  be  represented  by  the  potential 
4X(r,Q)  =  <)>oexp(ikdR)/R 

(4) 

■  -i  $  k .  £  (2n+l)h  (k.r  )j  (k  r)P  (cos0) 
odn  nas  nd  n 

where  k.  =  w/c.  and  r  and  r  refer  to  source  and  receiver  locations 
das 

respectively,  with  rg  >  1.  The  spherical  Hankel  and  Bessel  functions  are 

denoted  by  h  (x)  and  j  (x) ,  and  P  (x)  is  the  Legendre  function.  For 
n  n  n 

k^rg  »  i  the  large  argument  asymptotic  expression  for  the  spherical 
Hankel  function  may  be  used  and  Eq.  4  reduces  to  the  expansion  for  an 
incident  plane  wave 

^(r.S)  »  E  (2n+l)inj  (k.r)P  (cos0)  (5) 

o  n  n  d  n 

with 

♦o  "  *oexp(ikdrs)/rs 

Alternatively,  a  distortional  (axisymmetric)  source  may  be  specified. 


i  2 

X  (r,0)  ■  Xokgrssin0exp(ikgR)/R 


-ix  k  L  ( 2n+l) h  (k  r  ) j  (k  r) dP  ( cos0) / d0 
ssn  nss  ns  n 


where  kg  *  u/cg.  For  this  source  the  vibration  is  transverse,  and  in  the 
incident  plane  at  large  distances. 


B.  Formulations  for  Solution 
1.  Wave  Harmonics 

Using  spherical  wave  harmonics,  the  potentials  which  define 
the  scattered  field  for  the  (axisymmetric)  problems  posed  above  are  given 
in  series  form 
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<P  “  I  A  h  (k.r)P  (cos6) 
n  n  n  a  n 

X-IBh  (k  r)dP  (cos0)/d6 
n  n  n  s  n 


(7) 


The  unknown  coefficients  A^  and  are  determined  from  the  appropriate 

boundary  conditions  and  source  configuration.  For  example,  in  the  case  of 

9 

a  dilatational  source  incident  on  a  spherical  cavity 


with 


alla22-a12a21 


and  where 

a1L  *  [n(n-l)-(kga)2/2]  hn(kda)+2(kda)hn+1(kda) 

a.,  -  n(n+l)  [(n*l)h  (k  a)-()t  a|h  .  (k  a)  ] 

12  n  s  s  n+l  s 

*21  *  ln-1)h„<Vl'<kd*)hn+l(kda) 

*22  '  I"2-l-|ks*>2'21VV)*(V,hn+llks*> 


(10) 


b1  =  (n(n-l)-(kga)2/2]j^(kda)+2(kda) jn+1(kda) 

b2  "  (n-1)3n(kda)"(kda)jn+l(V) 

Substituting  Eqs.  10  into  Eq.  9 

D  =  Wkda)  Vksa)+Wkda)Vl(ksa) 

+*3hn+l (kda)  hn(ksa)  +^4hn+l  (kda)  hn+l  (ksa) 


(ID 


with 
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♦,  =  (k  a)2/2{(2n+l) (n-l)-(k  a) 2/2  -  4n(n-l) |n2-l)/(k  a)2} 
1  s  s  s 

♦2  =  (kga) [n(n-l) (n+2)-(kga)2/2] 

*3  =  (kda) [(n2-l) (n+2)-(kga)2] 

<j>.  =  -(k.a) (k  a) [ (n-1) (n+2) ] 

4  as 

Using  the  recursion  relation 

h  ,  (x)  =  (x)  -  d[h  (x)  2]/dx 

n+1  x  n  n 

The  quantity  D  can  also  be  expressed  in  the  following  form 

D  =  Lh  (k,a)h  (k  a)+I  h  (k.a)h’(k  a) 

In  d  ns  2nd  ns 

+* ,h ' (k  a) h  (k  a) +i .h‘ (k  .a) h’ (k  a) 

3  n  d  ns  4  n  d  ns 

with 

♦.  *  n(n2-l) (n+2)+(k  a)4/4  -  (k  a) 2 [2n(n+l)-l]/2 
1  s  s 

$2  =  (ksa) [(ksa)2/2] 

-  (kga) [ (ksa) 2- (n+2) (n-1) ] 
i4  *  (kga) (k^a) [2-n(n+l) 1 


(12) 


(13) 


(14) 


(15) 


2.  Creeping  Waves 

In  general  and  including  the  results  above,  the  wave  harmonic 
formulation  for  elastic  scattering  can  be  shown  to  yield  convergent  series 
for  the  potentials  as  well  as  for  the  associated  displacement  and  stress 
fields.  However,  it  is  the  nature  of  the  formulation  that  the  number  of 
terms  required  for  convergence  increase's  with  frequency.  In  order  to 
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transform  such  series  into  ones  which  converge  rapidly  at  high  frequencies, 
the  Watson  transformation  is  employed.  A  key  ingredient  in  this  trans¬ 
formation  is  the  formula 


£f (n+1/2)  =  (1/2) 
n 


/ 


f  ( v)  exp  (-iirv)  /cosirvdv 


(16) 


where  f  is  a  reasonably  well  behaved  function  and  C  represents  an 
appropriate  contour.  This  contour  is  then  deformed  such  that  the  integral 
can  be  evaluated  in  terms  of  a  set  of  pole  contributions,  each  representing 
a  creeping  wave.  Specifically,  for  the  class  of  problems  considered  here, 
this  approach  leads  to  the  following  identity 


Zn ( 2n+l) P  (cos0) f (n) /D(n) 
n=u  n 

=*  2tt^{vPv  1^2(cos(tt-0)  )  f  (v-1/2)  /  fcosvir3D(v-l/2)/3v]  }v=v 

where  the  (creeping  wave)  roots  satisfy  the  equation 
D( Vj-1/2)  =  0 

(It  should  be  noted  that  Eq.  17  may  also  be  obtained  by  means  of  the 
Poisson  summation  formula.11-) 


(17) 


(18) 


For  example,  we  again  consider  the  case  of  a  dilatational  source 

incident  on  a  spherical  cavity.  The  solutions  to  Eq.  18,  for  k  a  »  1,  are 
12  13 

of  three  types.  '  An  infinite  set  of  shear  type 


,  .  _  -1/3  5/3iri,  -iir  ,  .1/2  ,  „  , 

v. .(n  )~n  +2  e  (e  »  3  =  1,  2,  3, 

X  J  S  S  J  s 


(19) 


An  infinite  set  of  dilatational  type 

±  v_.(fl.)'fl.+21/3e5/3,ti(eiV>  (fi.)1/3  ,  j  =  1,  2,  3,  ...,  (20) 

2j  a  a  3  a 

and  a  single  pair  of  the  Rayleigh  type 
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where  the  imaginary  part  of  v  vanishes  exponentially  as  (}-+“.  The  z. 

R  14  b  3 

are  the  roots  of  Ai(z)=0  which  are  real  and  negative,  Ai(z)  is  the  Airy 

function  and  c  is  the  Rayleigh  wave  speed.  It  is  interesting,  as  well 

K 

as  reassuring,  to  see  that  for  high  frequency  scattering  from  a  cavity 

there  is  a  root  v  -  uja/c  which  corresponds  to  the  Rayleigh  surface 
R  R 

wave  on  a  (planar)  half  space. 

13  15 

The  corresponding  potentials  (for  0^0)  are  given  by  ' 


d>  =  2iri  exp(-iir/4)  $  I 
Tc.w.  e  os 


(-1)S  ?v  Ij_(kdr) 

3  3  v 


+  (-ni/2)exp(iirv.)  (a?2b1~a13b2) 

3  “  n=v 


(22) 


h  (k  r)  3/OD/3v) 

-  Q 

V 


P  (cos8)  exp  [iv  .  (2tts+7t/2)  ) 
n=v  v  3 


Xc  w  =  2iriexp(-iTr/4)  |  (-1)  (iri/2) expfiirv.. )  (a^bj-a^b^ 


n=v 


(23) 


i  (k  r)/(3D/3v) 

V j-1/2  3 


dP  (cos9)/d0 
n=v  v 


with  v  =  v_.-l/2  and  0^0.  The  travelling,  or  creeping,  wave  nature  of  the 
above  potentials  is  evident  if  one  substitutes  the  large  order  asymptotic 
expressions  for  the  Legendre  functions  in  Eqs.  22  and  23. 


P  (cos0)  -  (27TVjSin0)  1/,2{exp[i  (Vj0-rr/4)  ]+exp[-i(Vj9-u/4)  ]  } 


(24) 


It  can  be  shown  that  for  large  f 2.  and  0  the  above  series,  as  well  as 

Q  s 

-N 

those  for  the  corresponding  displacements,  converge  as  v.  with  N  >  1  for 

13  3 

0  <  0  <  tt/2  when  s=0,  and  for  0  <  0  <  it  when  s  f  1.  Thus,  unlike  the 
wave  harmonic  series  given  earlier  the  convergence  of  the  creeping  wave 
series  actually  improves  with  increasing  frequency  within  their  domain  of 
validity.  Also,  convergence,  improves  with  increasing  9.  In  other  words, 
it  improves  as  one  moves  away  from  the  illuminated  zone.  This  is  a  general 
phenomenon  which  in  some  cases  limits  the  utility  of  creeping  wave 
solutions  to  the  shadow  zone  therefore  providing  a  convenient  supplement  to 
physical  optics  solutions  (Section  II. B. 3). 


-8- 
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Creeping  wave  solutions  exist  for  inclusions  of  other  than  spherical 

geometry  and  zero  impedance  as  well  as  for  other  types  of  incident  waves 

both  harmonic  and  transient.  Cylindrical  cavities  have  been  studied  by 
16  17 

Miklowitz  and  Peck  .  Gaunaurd  and  Uberall  have  obtained  solutions  to 

many  different  scattering  problems  including  fluid  filled  spherical  cavities 

18  19 

embedded  in  a  viscoelastic  as  well  as  elastic  medium.  '  (In  their 

approach  the  contributions  of  creeping  waves  are  represented  in  terms  of 

Regge  poles.)  They  have  also  considered  incident  waves  of  the  distortional 
20 

type. 

3.  Physical  Optics;  Kirchhoff 

The  physical  optics,  or  Kirchhoff,  high  frequency  asymptotic 

21  22  23 

approach  to  scattering  has  been  used  extensively  for  radar  and  acoustic  ' 
applications.  The  method  is  applicable  to  the  illuminated  zone  and  is  most 
useful  when  applied  to  the  specular  return.  The  basis  for  this  formulation 
is  that  each  differential  scatterer  acts  as  though  it  were  located  in  a 
planar  infinite  baffle.  For  the  problem  of  elastic  scattering,  this 
assumption  allows  one  to  express  the  scattered  potentials  in  terms  of  a 
Green's  function  which  satisfies  Neumann  boundary  conditions  at  the 
boundary  of  the  scatterer. 


4>S(R)  =  -  2  /  G  (R, R  )  3d>2(R  )/3n  dS 

J  y  o  o  o 

s 

o 


XS(R)  =  -  2  /  G  (R,  R  )  3XS(R  )/3n  dS 

J  X  o  o  o 

s 

o 


(25) 


with  G  (R,R  )  =  exp(ik,  |R-R  |)/[4tt|r-R  I] 
<p  o  d 1  o  o' 

Gx(R,Rq)  =  exp(iks|R-Ro|)/[4Tr|R-Ro|) 


(26) 


and  where  Sq  is  the  surface  of  the  scatterer,  n  represents  the  normal  to 
the  surface  and  R  denotes  a  vector  quantity.  The  integrals  in  Eqs.  25 
when  evaluated  asymptotically  for  large  k  a  (or  k  a)  by  means  of  stationary 

Q  S 
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phase,  yield  the  specularly  scattered  field  for  either  bistatic  or  mono- 
static  conditions.  The  values  for  the  normal  derivatives  of  the  potentials 
evaluated  at  the  surface  are  obtained  from  reflection  coefficients 
associated  with  the  corresponding  planar  scattering  problem.  For  example, 
for  the  problem  at  hand,  consider  a  plane  dilatational  wave  incident  on  a 
free  surface  at  angle  0^  from  the  normal. 

*  fcexptik^fxsinS^-zcosB^)  ]  (27) 

24 

Both  a  dilatational  and  a  distortional  wave  are  scattered 

<)>s  *  A$exp[ik,(xsin0 .  +  zcos0.)] 
d  l  l 

s 

X  =  B$exp [ikg (xsin6  +  zcos5) ] 

with  k  ,sin0 .  =  k  sin6 
d  l  s 

and 

2 

c  cosO.tan  2<5  -  c.cos6 
si  d 

- - 2 - 

c  cosG.tan  26  +  c.cos6 
si  d 

c  sin20.  2c  cos  6 

_  _s _ i  _ s _ 

c.cos2  6  .  ^  2„  .  . 

d  cos 0. tan  2  6  +  c.cos6 

l  d 

Therefore,  from  Eqs.  28 

d<t>S/dn  =  d(fiS/dz  =  ik.COS0,4>S 

d  1 

and  (30) 


(28) 


(29) 


s  s 

3x  /3n  ■  ikgcos6x 

Thus,  under  Kirchhoff  approximations,  although  mode  conversion  is  accounted 
for  (and  neglecting  multiple  scattering) ,  the  scattered  potentials  can  be 
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~  ’'T' . 


considered  independently.  Actually,  if  one  limits  oneself  to  back 
scattering,  then  symmetry  considerations  dictate  that  the  exact  solution 
as  well  yields  no  coupling. ^ 

For  the  problem  of  a  plane  dilatational  wave  incident  on  a  spherical 
cavity  the  dilatational  component  of  the  (bistatic)  scattered  field  is  given  by* 

g  -ik^acosy  ikrfp 

<t>K  =  •facosysinye  e  /5  (31) 

2  2  1/2 
where  p  =  [F  a  -2aRcos (y-0) ] 

3  1/2 

C  =  pRsinysinQ [cosy+(a/R) sin  y/sin(y-0) -siny/tan(y-0) ]  ' 

and 

(R/a) sin(2y-0)  =  siny 

In  the  far-field,  i.e.,  for  a/R  <<  1,  Eq.  31  reduces  to  the  well  known 
result  for  specular  scattering  (both  bistatic  and  monostatic) 

-  -(a/2Rleik<r-acose/2>  (32) 


or 


R<t>^/  4>i 


=  a/2R 


(33) 


It  is  the  nature  of  this  formulation  that  it  is  easily  constructed  for 
arbitrary  geometry  and  (locally  reacting)  surface  impedance. 

4.  T-Matrix 

The  T-matrix  formulation  of  scattering  has  significant 
computational  advantages  over  other,  more  traditional,  matrix  formulations. 
These  advantages  relate  to  the  nature  of  the  elements  of  the  matrix  which 
turn  out  to  be  surface  integrals  exhibiting  regular  integrands.3  The  first 
step  in  the  formulation  is  to  express  the  incident  and  scattered  fields  in 


‘This  result  was  obtained  by  my  colleague  Dr.  W.  T.  Ellison  and  represents 
an  extension  of  his  previous  work  on  monostatic  scattering. 26 
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terms  of  the  appropriate  eigenfunction  series  for  incoming  and  outgoing 

waves,  viz.  cylindrical  harmonics  for  two  dimensional  scattering  problems 

and  spherical  harmonics  for  three  dimensional  problems.  Then,  a  set  of 

functions,  f^,  is  introduced  which  is  complete  but  may  or  may  not  be 

orthogonal,  on  the  scattering  surface.  This  set  of  functions  is  not 

necessarily  unique  and  may  be  chosen  for  analytical  or  computational 
4 

considerations.  Finally,  the  incident  and  scattered  fields  are  expanded 
in  terms  of  these  functions  f ^  on  the  surface  of  the  boundary  and  the 
appropriate  boundary  conditions  are  satisfied  term  by  term.  The  T-matrix, 
which  relates  the  unknown  coefficients  of  the  series  representing  the 
scattered  field  in  terms  of  the  known  incident  field,  is  thus  derived. 

For  completeness  the  procedure  for  scalar  scattering  is  outlined  below. 

Consider  three-dimensional  scalar  scattering  and  a  scatterer  that 
satisfies  Dirichlet  boundary  conditions.  The  incident  wave  is  given  by 


<fi  (r,d,ip)  =  g  AnCn(r,6,4>) 


(34) 


The  scattered  wave  is  represented  by 


<)>  (r , 0, 4>) 


E  B„n_(r,0,i)i) 
m  mm 


(35) 


For  example,  for  axisymmetric  problems 
f,n(r,0,$)  =  jn(kr)Pn(cos0) 

and 

n  (r ,0,<p)  *  h  (k  )P  (cos0) 
n  n  r  n 

where  j  and  h  denote  the  spherical  Bessel  and  Kankel  functions  respectively 
n  n 

and  Pn  is  the  Legendre  function.  The  Dirichlet  condition  on  the  boundary  Sq 


is  given  by 


(♦*  + 


(36) 
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Expanding  Eqa.  34  and  35  in  teraa  of  the  functions  f ^  and  satisfying  Eq.  36 
yields  the  set  of  equations 


i'VW  *  W'"” 


where  (u,v) 


/ 

S 

o 


u*vdS 

o 


n* j  -  0,  1,  ...» 
(r,6,t)  on  Sq 


(37) 


(37a) 


Using  the  usual  notation,  Eqs.  37  (truncated)  may  be  written  in  matrix 
form 

[Q]A  +  [Q]S  -  0  (38) 


where  A  * 

[A^  $  •  “ 

"V 

5  - 

•”bn] 

[Q] 

'  «jn  ' 

tQl 

"  QJn  “ 

(frnn) 

Solving  Eq.  38 

5  -  - [Q]'1 [Q)A  *  TA  (39) 

where  T  defines  the  T-matrix  which  in  turn  defines  the  scattered  field  as 
given  by  Eq.  35. 
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j 

i 

III.  THE  USE  OF  HIGH  FREQUENCY  ASYMPTOTICS  TO  EVALUATE  AND  FACILITATE  THE  j 

CONVERGENCE  OF  T- MATRIX  SOLUTIONS  ; 

A.  General  Formulation  j 

As  indicated  above,  in  many  cases  high  frequency  scattering  solutions 

j 

share  the  simplicity  often  found  in  low  frequency  formulations.  For  example,  j 

while  in  the  low  frequency,  or  Rayleigh,  region  only  the  volume,  density 

and  compressibility  of  a  complex  scatterer  may  be  needed,  in  the  high  i 

j 

frequency  regime,  often  only  the  local  radii  of  curvature  and  reflection 

coefficient  need  be  known  to  calculate  the  specular  return.  Thus,  as  is  | 

.1 

the  case  for  low  frequencies  the  nature  of  the  solution  to  a  variety  of  j 

complex  scattering  problems  tends  to  coalesce  yielding  a  large  number  of 
scattering  problems  which  is  bounded  at  both  the  low  and  high  frequency 
ends  by  known  solutions.  These  solutions  can  benefit  T-matrix  formulations 
in  a  number  of  ways.  First,  they  can  provide  convergence  checks.  These 
checks  need  not  necessarily  be  in  terms  of  matching  amplitudes  but  also  in 
the  case  of  creeping  wave  solutions,  the  checks  may  be  in  terms  of  phase 
velocities  (Eqs.  22  and  23)  or,  for  transient  problems  arrival  times. 

Also,  although  only  separable  geometries  have  been  considered  in 
Section  II,  the  asymptotic  results  can  be  extrapolated  to  other  more 
complex  geometries.  This  is  straightforward  for  the  specular  return  since 
it  is  a  local  phenomenon.  However,  as  stated  above,  it  is  also  the  case, 
albeit  somewhat  more  limited,  for  the  creeping  wave  results  by  proper  inter¬ 
pretation  of  the  functional  dependence  of  the  creeping  waves  on  the  spatial 
coordinates . 

Additionally,  another  use  for  these  solutions  is  being  proposed 
here.  In  particular,  it  is  suggested  that  there  may  be  computational 
advantages  to  be  gained  by  integrating  these  solutions  into  the  T-matrix 
formulation  itself. 

One  possibility  would  be  to  use  the  high  frequency  solutions  as  an 
aide  in  choosing  appropriate  base  functions,  f y  As  mentioned  earlier, 
the  base  functions,  ty  can  be  chosen  for  analytical  considerations, 
e.g.,  consistency  with  the  Helmholtz  integral  formulation,  as  well  as  for 
computational  considerations,  e.g.,  minimization  of  the  mean  square 
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deviation  from  the  prescribed  boundary  conditions.^  In  terms  of  this 
latter  consideration,  it  would  be  reasonable  to  expect  that  knowledge  of 
the  nature  of  the  high  frequency  behavior  at  the  boundary  should  allow 
one  to  develop  base  functions,  either  in  part  or  the  complete  set,  which 
would  yield  computational  benefits  at  high  frequencies.  Unfortunately, 
this  approach  would  generally  lead  to  frequency  dependent  base  functions 
and  it  is  not  clear  to  the  author  how  this  could  be  conveniently  accompli shod 
in  a  general  fashion.  Thus,  this  possibility  is  not  pursued  further. 

However,  what  is  suggested  for  further  study  is  the  following. 

Consider  the  T-matrix  formulation  as  presented  in  the  previous 
section.  Let  represent  the  (known)  high  frequency  asymptotic  solution 
for  the  scattered  field  which,  like  the  actual  field,  can  be  expanded  in 
terms  of  two  or  three  dimensional  harmonics  for  outgoing  waves. 

■  SW-9-*1 

with 

c»  “ 

Thus,  since  nm  are  complete  over  all  space  V,  Eq.  35  may  be  written  as 

$S(r,e,$)  -  ♦^(r,e,$)  ♦  ♦  (r,e,*)  (42) 

with 

•L11'6'*1  <«> 

The  matrix  equation  (Eq.  37)  becomes 

S<\  <f3-Cn)  +Dn(fj,nn)J  -  *  (r,e,*)  on  So  (44) 

where  0  =  B  -C 

n  n  n 


In  matrix  notation 
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with  i*  - 

What  is  accomplished  by  the  above  manipulations  is  that  the  (unavoidable) 
computational  difficulties  that  are  associated  with  the  series  expansion 
of  the  scattered  field  at  high  frequency  have  been  isolated  and  extracted 
from  the  T-matrix.  (It  should  be  noted  that  this  approach  is  equally 
applicable  to  the  wave  harmonic  formulation  for  separable  geometries . ) 

B.  Illustrative  Examples 

1.  Spherical  Cavity 

To  illustrate  the  concepts  discussed  in  Section  A  we  first 
continue  our  consideration  of  scattering  from  a  spherical  cavity  and  in 
particular,  backscattering.  For  simplicity,  we  will  limit  our  high 
frequency  asymptotic  solution  to  the  specular  return  given  in  Section 
II. B. 3.  In  other  words,  let  be  given  by  Eq.  31.  Further,  since  for  a 
spherical  scatterer,  nothing  is  gained  by  choosing  the  set  f ^  to  be  other 
than  the  Legendre  functions  P^ (coa6) .  Results  for  the  backscattered 
pressure  and  relative  phase  vs.  frequency  are  shorn  in  Figs.  2  and  3.  The 
exact  solution  was  obtained  by  summing  the  wave  harmonic  series  (Eqs.  7- 
10),  or  equivalently  the  T-matrix  series  (Eqs.  35  and  37),  to  convergence. 
Roughly  speaking,  for  the  problem  at  hand  convergence  requires  N  >  2(ka+l) 
The  dashed  line  connects  points  which  represent  the  contribution  of  the 
first  term  in  the  T-matrix  series  which  was  calculated  using  only  a  2  x  2 
matrix.  Obviously  the  result  is  adequate  only  for  low  values  of  ka.  The 
dotted  line  is  the  corresponding  result  for  the  modified  T-matrix 
formulation  (Eq.  42) .  Further  improvement  could  be  obtained  by  including 
the  creeping  wave  (Eq.  22),  as  well  as  the  specular,  contribution  to 

2.  Prolate  Spheriodal  Cavity 

The  example  of  a  spherical  scatterer  is  of  sosMwhat  limited 
value  as  an  illustrative  example  in  the  sense  that  the  T-matrix  approach 
degenerates  identically  to  the  wave  harmonic  solution.  Thus,  in  this 
section  we  consider  (scalar)  scattering  from  a  cavity  in  the  shape  of  a 


(46) 
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prolate  spheroid.  For  convenience,  we  again  limit  ourselves  to  back- 
scattering  and  consider  end-on  incidence  (Fig.  4) .  The  prolate  spheroidal 
coordinates  (£,n)  are  related  to  cartesian  coordinates  by  the  transformation 

x  -  (d/2)  (U-n2)  (C2-i)l1/2 

(47) 

z  -  (d/2)n£ 

with  -1  <  n  <  1  and  1  s  {  <  ■ 

The  scattering  surface  is  defined  by  the  equation  £  -  £  .  For  the  present 
problesi  £q  is  taken  to  be  1.1547,  which  yields  a  spheroid  with  a  2:1  aspect 
ratio.  For  this  geometry,  and  assuming  symmetry  about  z,  the  inner  product 
<£q.  37a)  becomes 

1 

(u,v)  -  2x(d/2)2(£2-l)1/2  f  (u*r) (£2-n2)1/2dn  (48) 

o  J  o 

-1 

Also,  the  local  radius  of  curvature  at  (£o»l),  i.e.,  end-on,  is  given  by 

a*  -  (d/2) (£2-l)/£  (49) 

o  o 

Thus,  the  high  frequency  (Kirchhoff)  backscattered  level  (equivalent  to 
Eq.  33  for  the  spherical  scatterer)  is  given  by 

|A4*/4i|  -  a*/2R  (50) 

29 

A  comparison  between  this  asymptotic  result  and  the  exact  solution  is 
shown  in  Fig.  5.  In  order  to  maintain  a  notation  consistent  with  Reference 
29  the  ordinate  in  Fig.  5  is  given  in  terms  of  a  backscattered  cross-section 
which  is  defined  by 

a/ak  -  l2R/(d/2)n£o/(£2-l))[^i)  2  (51) 

Also  shown  in  Fig.  5  is  the  result  of  a  T-matrix  formulation  using  base 
functions 

f^  -  cost je)  j  -  0,  1,  ...  N  (52) 

where 
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-1  -125  1/2 

0  =■  tan  (x/z)  -  tan  {[(1-n  )(E  -Dl  /(nE  )  (53) 

o  o 

Consistent  with  the  previous  example  the  number  of  terms  required  for 
convergence  is  given  by 

N  =  2 (k£Qd/2+l)  (54) 

The  integration  implied  by  Eq.  48  was  performed  using  the  trapezoid  rule 
with  a  mesh  size  equivalent  to  A/8. 

It  should  be  noted  that,  as  was  the  case  in  the  previous  section, 
the  high  frequency  asymptotic  solution  is  reached  before  the  convergence 
criterion  (Eq.  54)  becomes  prohibitive.  An  indication  of  the  effectiveness 
of  the  high  frequency  modification  for  this  case  is  shown  in  Fig.  6.  Here 
the  percent  error  incurred  by  constructing  the  scattered  field  using  only 
the  first  N/2  terms  in  the  T-matrix  solution  is  plotted  vs  non-dimensional 
frequency  for  the  modified  as  well  as  the  unmodified  formulation. 

In  summary,  both  examples,  although  elementary,  indicate  that  the 
extraction  of  appropriate  high  frequency  components  from  acoustic  (or 
elastic)  scattering  problems  may  represent  a  viable  means  for  improving 
the  convergence  of  T-matrix  formulations  to  more  complex  problems. 
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[h  Frequency  Modification  to  T-Matrix  Formulation. 


Relative  Phase  of  Backscattered  Potential  from  a  Spherical  Cavity  vs.  ka 
Effect  of  High  Frequency  Modification  to  T-Matrix  Formulation. 


Geometry  for  Scattering  from  Prolate  Spheriod. 
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fig.  5  -  Magnitude  of  backscattered  cross-section  from  a  prolate  spheriodal 
cavity  vs  kf,  d/2. 
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Percent  error  in  backscattered  cross-section  for  a  prolate  spheroidal  cavity  using 
terms  for  $  (N  given  by  Eq.  54) . 


